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Abstract 

We analyze a simple stochastic model to describe motor molecules which 
cooperate in large groups and present a physical mechanism which can lead to 
oscillatory motion if the motors are elastically coupled to their environment. 
Beyond a critical fuel concentration, the non-moving state of the system be- 
comes unstable with respect to a mode with angular frequency uj. We present 
a perturbative description of the system near the instability and demonstrate 
that oscillation frequencies are determined by the typical timescales of the 
motors. 

PACS Numbers: 87.10.+e, 05.40.+j 
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Motor proteins are highly specialized macromolecules which can consume mechanical 
energy to induce motion and to generate forces. These molecules are involved in active 
transport processes, cell locomotion and muscle contraction . A typical motor molecule 
specifically attaches to a certain protein filament which serves as a track for its motion. In 
the presence of fuel, which in the cell ist adenosinetriphosphate (ATP), the motor starts 
moving in a direction defined by the polarity of the track filament. Experimental methods 
to study the physical properties of motor proteins have been developed in recent years ||-f7| . 
These methods allow to measure forces and velocities of individual motors or small groups 
of motors. 

Normally, at given fuel concentration and temperature, motor molecules will generate a 
constant average force leading to a constant average velocity JFj. In some cases, however, 
biological motors are used to drive oscillatory motion. Insect flight muscles e.g. generate 
oscillating forces in order to move the wings with high frequency. While for some insects (e.g. 
butterflies) the timing of these oscillations is controlled by a periodic external nerve signal, 
others (e.g. bees and wasps) generate oscillations within the muscle. The mechanisms which 
lead to these oscillations are not understood but they seem to be related to the process of 
force-generation ||. 

The purpose of this article is to describe theoretically a physical mechanism by which 
motor collections acting on a spring can give rise to oscillatory motion. Oscillations occur 
as the result of cooperation of many motors while individual motors only generate time- 
independent velocities. 

The microscopic mechanisms which are responsible for the force generation of molecular 
motors are not understood. In order to indentify possible physical mechanisms which allow 
such a motor molecule to convert the chemical energy of a fuel to motion and mechanical 
work, simple physical models have been suggested [|9|-|T6||. One idea |TIJ is to simplify the 
internal degrees of freedom of the motor to two different states of a particles which moves 
along a one-dimensional coordinate x. The interaction of the particle with the track is 
described by periodic potentials which reflect the periodic structure along the surface of the 
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track. Within such a two-state model motion is induced if the potentials are asymmetric 
with respect to x —>■ —x (which reflects the polarity of the track) and if detailed balance of 
the transitions between states is broken 111 Oil. 



New phenomena occur if the collective motion of many particles is considered |TT|-|T9 
In a simple model for collective motion, particles are rigidly coupled to a common structure 
which they set in motion collectively |18| . As a result of cooperation, motion can in this 



case even occur in a symmetric system via spontaneous symmetry breaking. Instabilities and 
dynamical phase transitions also occur in asymmetric systems revealing the rich behaviour of 
many-motor systems. At these transitions, a steady state becomes unstable and the system 
chooses a new steady state. These transitions correspond therefore to dynamical instabilities 
at angular frequency uj = 0. 

In the following, we demonstrate that many coupled motors also can induce instabilities 
with uj ^ towards a state with periodically varying velocity. Such instabilities occur if 
the system is elastically coupled to its environment. This implies, that the complete system 
is connected to its environment via a spring attached to the common backbone, see Fig. [I] 
|20|| . A single particle reaches in this case a steady state with zero velocity and a position 
at which the average force generated by the motor is balanced by the elastic force of the 
spring p. 

We are interested in the behavior of many motors coupled together. Let us first review 
the most simple two-state model which we are going to discuss fllOyiSf . Motors are described 
by particles which move in two periodic potentials W a (x) with period I. Transitions between 
states a = 1,2 occur with rates uji{x) and u^O^) which also are /-periodic. The excitation 
amplitude can in general be written as [[U]] u>i(x) = u>2(x)[exp[(Wi — W2)/T] +Q9(x)], where 
T denotes temperature measured in units of the Boltzmann constant. We assume that the 
periodic function 9 (x) is given. In principle it can be calculated from a specific model for the 
reaction kinetics of the fuel coupled to the particle []2T|. The amplitude Q is related to the 



difference of chemical potentials A/i = hatp — Hadp — Hp of the fuel and its products which 
is the chemical driving force of the motors. Here, we assume that the fuel is ATP which 



is hydrolyzed to adenosinediphosphate (ADP) and phosphate (P): ATP ADP + P. For 
A/x = 0, Q = and the rates u)\ and u 2 obey detailed balance uji/u 2 = exp[{W\ — W 2 )/T]. As 
soon as A/x 7^ 0, f2 becomes nonzero. For small A/x/T, Q ~ A/x; for large A/x/T, f2 ~ Catp 
is proportional to the fuel concentration. 

We now study the behavior of many particles which are rigidly connected to a common 
backbone, see Fig. |l[ Particles have a fixed spacing q which for simplicity is assumed to be 
incommensurate with the period I, i.e. l/q is irrational. In the limit of an infinite system 
one can then introduce densities P CT (£) with a = 1, 2 which give the probability to find a 
particle at position £ = a; mod / relative to the beginning of the potential period in state a. 
These densities are not independent but obey Pi(£) + P%{£) = l/l and J" d£(Pi + P 2 ) = 1- 
The equations of motion for this system read Jl8| 

d t P x + vd^P x = -W1P1 + cu 2 P 2 

a t P 2 + W^P 2 = WiPx - C^ 2 P 2 (1) 

f ext = Xv + KX + f diPxd^Wx - W 2 ) . (2) 

Eq. (|I]) describes the dynamics of the density Pi(£) resultig from motion of the backbone 
with velocity v = dtX and the transitions between the states. Eq. (0) is a force balance: 
The externally applied force per particle f ex t{t) is balanced by viscous drag with damping 
coefficient A, the average force excerted by the potentials and an additional elastic force KX 
of a spring of length X with elastic modulus K N where N is the number of particles. Note, 
that Eqns. (|l|) and (§) only depend on the difference W\ — W 2 of the two potentials and we 
can choose W 2 to be constant without loss of generality. 

As a consequence of the spring action, a nonmoving solution to Eqns. (0) and (0) 
exists for f ext — 0, v — 0, Pi = R = u 2 /(al), where = ct>i(£) + cu 2 (£) and X = 

X = — /q d£i25^(Wi — W%)/K. We first study the linear stability of this solution and 
determine the instability threshold with respect to oscillations. With the ansatz Pl(£, t) = 
R(£) + p(0 exp(st), v(t) = uexp(st) and X(t) = X + wexp(st)/s, one finds using Eq. ([I]) 
and @ to linear order in p and the velocity amplitude u: 
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P(0 = -u 



(3) 



s + a(0 



The possible values of the complex eigenvalue s = — r + iu are determined by 




•* dzRdz(W 1 -W 2 ) 
i s + a(£) 



(4) 



The non-moving state is unstable if r < 0. The instability occurs for r = 0, where the real 
part of s vanishes. This happens at a threshold value Q = VL C (K) for which a(£) = a c (0- 
This critical value together with the frequency uo c at the instability is determined by 



In the limit K = where no elastic element is present, it follows from Eq. (§) that uo c = 
0. Eq. (0) now determines Q c (0) which is the condition for the instability of the non- 



instability leads to spontaneous motion via a symmetry-breaking transition. For nonzero K 
the instability occurs for £l c (K) = f2 c (0) + SQ C with 5Q C ~ K for small K and the system 
starts to oscillate with finite angular frequency uj c ~ K 1 ! 2 . Note, that beyond a maximal 
value K > K max , the resting state is stable for any value of Q. 

Fig. |2| (a) and (b) show examples for these oscillations which have been calculated numer- 
ically for constant deexcitation rate u> 2 , \uJ2l 2 /U = 0.1 and piecewise linear and symmetric 
potentials W a (£), see Fig. [I] (a) with a = 1/2. In this case K max l 2 /U ~ 2.8. The function 
is chosen as shown in Fig. [I] (a) with d/l = 0.1. Fig. |2] (a) displays the position X 
versus time t for an excitation level Q = 0.1 slightly above threshold and elastic modulus 
Kl 2 /U = 0.2. After an initial relaxation period, motion is almost sinusoidal. An example 
for small elastic molulus Kl 2 /U = 0.002 and same Q is shown in Fig. (b). The periodic 
motion in this case shows cusp-like maxima which correspond to sudden changes of the 
velocity. These cusps are related to discontinuities of the velocity as a function of external 




(5) 



(6) 



moving state as previously described in Ref. ||18|| . In the case of a symmetric system, this 



force for K = || . 
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We now give an analytical description of oscillations in the vicinity of the instabil- 
ity. Anticipating, that the motion is periodic with period tp = 2tt/uj, we can write 
Pi(t,t) = ET=-ooPi(^k)e^\ v(t) = Tk=-ooVke ik "\ and f ext (t) = EZ-oo fke lhwt , which 
defines the Fourier coefficients Pi(£, k), Vk and Using this representation, one can derive 
the nonlinear relation between velocity and external force 

/* = F M V l + F MmVlV m + F^VlVmVn + 0(V) . (7) 

The coefficients F^j _ kn can be calculated by first rewriting Eq. ([[]) as 

Pl (^k) = S k>0 R(0-j:-^TV l d^,m) . (8) 

Inserting the ansatz 

P 1 (^k) = R5k, + P^(0vi + PS ) m(0viVm + O(v 3 ) , (9) 

into Eq. (H), one obtains a recursion relation for the functions P k k 1 ..k n ' 

p«> = - E ^^fcl, (io) 



, a + iuok 



Using Eq. (0), one finds 



y iujk Jo a + iuk J 

The linear response function of the system is given by the inverse of F k j\ Therefore the 
instability condition Eq. (§) corresponds to F^j' = 0. For n > 1, 

Fit,..,K = f d£P (n \OkM,,kA(Wi ~ W 2 ) . (12) 

<j 

Note, that the coefficients kn are nonzero only if A; = ki+..+k n . This results from time 



translational symmetry [22]. If no elastic coupling to the environment is present, i.e. K — 0, 
one recovers for constant external force fo the steady states which have been described 
previously [0: f ext = fo(vo); v n = for n ^ 0. As soon as K ^ 0, vq must vanish and a 
constant external force /q only changes the average position Xq. 



Spontaneous oscillations are solutions to Eq. (|7|) for /„ = 0. The dominant terms near 
the instability of V\ are given by 

= Fffvi + G®v- X V2 + G (3) ^_! (13) 
= F£ ) v 2 + Fg ) 1 v 2 1 , (14) 

where G (2) = f[%_ x + f[ 2) _^ 2 and G (3) = F^l x _ x + i^t-i,! + Fj® w . As soon as v x and v 2 
are determined, higher orders v n can be obtained recursively using Eq, (|7|). From Eq. (|i~4"D 
one finds that t>2 ~ f 2 . Inserting this value in Eq. (fL3f), one obtains 

= F^V + G^ 2 ^ , (15) 

with an effective third order coefficient 

£0) = G 0) _ F^G^/F^. One solution to Eq. 
(|I~5|) is always v-y = 0. The remaining solutions are described by 

H 2 = -F^/G^ . (16) 

Since the amplitude \vi\ 2 is a real number, solutions to Eq. (|H|) exist only if the right hand 
side of Eq. (|T6"D is real and positive. In general, however, F$ and G^ are complex numbers. 



The requirement to find real solutions to Eq. flTED fixes the oscillation frequency. This can be 
seen as follows: If the system is at its instability threshold Q = Q C (K), Fil = for u = u c , 
see Eq. (|J). In the (Q, a;)-plane, a line exists along which Fq /G^ is real. This line can be 
parametrized by a function uj = u s (Q) which at the instability threshold obeys uo c = u s (Q c ). 
For > tt c , —Fy^/G^ is real and positive along this line thus allowing for a solution to 
Eq. (|T6"D. For Q < Q c , v\ = is the only solution. Fig. ^| (c) shows the dependence of the 
selected frequency u s as a function of the excitation amplitude. For Q < Q C (K), u s = 0. At 
the instability, a finite frequency is selected which increases for increasing Q. 

Note, that this example is developed for a symmetric system for which an additional 
restriction holds: Eq. ([7|) is now invariant with respect to the transformation v n —>■ —v n , 
f n — > —f n . As a result, all even expansion coefficients in Eq. (g) vanish, F^ = = .. = 0. 
In particular, G^ vanishes and G^ = G^ 3 \ Furthermore, all terms which couple v± to modes 



v n with even n vanish and v-i = V4 = .. = 0. As a result, the periodic motion of a symmetric 
system has the symmetry property v(t + tp/2) = —v(t). 

The parameters used in Fig. |2] correspond e.g. to the choice A ~ 10~ 8 kg/s, UJ2 — 10 3 s _1 , 
/ ~ 10 _8 m, U ~ lOksT. Oscillation frequencies in this case vary between u s = for K — 
and uo s ~ buj2 for K ~ K max ~ 10 _3 A^/m per motor |]23[| . This shows that using typical 
time scales of biological motors the mechanism described here generates frequencies up to 
the kHz range. 

In summary, we have shown using a simple two-state model for foce-generation of molecu- 
lar motors that oscillatory motion can occur as a collective effect. Motors which individually 
only generate constant forces, show oscillating behavior if they are coupled in large numbers 
and work against an elastic spring. Beyond a critical amplitude of excitations which drive 
the system, a nonmoving state becomes unstable with respect to periodic motion. This ar- 
gument suggests that oscillations are a possible mode of motion of biological motor proteins 
and gives an example of how some insect flight muscles could generate oscillations via the 
mechanism of force generation. 
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FIGURES 

FIG. 1. (a) Periodic potentials W\ and W2 used to calculate the behavior of the system. The 
Function 9{x) is chosen to be nonzero in the vicinity of the potential minimum, (b) Many particles 
coupled rigidly to a backbone which is connected with its environment via a spring K. The particle 
spacing is denoted q. 

FIG. 2. (a) Position X versus time t for a symmetric system with d/l = 0.1, Vt = 0.1, 
\uJ2l 2 /U = 0.1 and Kl 2 /U = 0.2 (see Fig. [l] (a)), (b) Position versus time for same system 
but Kl 2 /U = 0.002. (c) Selected frequency uj s as a function of the excitation amplitude 17 for the 
case Q = 0.1, Kl 2 /U = 0.2. The critical excitation amplitude is f2 c ~ 0.08. 



11 




Fig. 1 (Julicher & Prost) 
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Fig. 2 (a) (Julicher & Prost) 




Fig. 2 (b) (Julicher & Prost) 




Fig. 2 (c) (Julicher & Prost) 



